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AN  OPTIMUM  DEPLOYMENT  OF  DEFENSIVE  WEAPONS 


ABSTRACT 

Dynamic  programming  methods  are  used  to  determine  the  allocation  of 
defensive  weapons  to  target  areas  so  as  to  minimize  effects  of  optimal 
attack  policies  of  the  offense. 
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INTRODUCTION 


The  introduction  of  long  range  ballistic  missiles,  with  their  large  de¬ 
structive  capabilities,  in  the  arsenals  of  war  weapons  introduces  some  new 
problems  in  planning  the  defense  of  key  industrial  and  military  sites.  The 
attack  times  involved  now  are  such  that  the  defense  weapons  cannot  be  rede¬ 
ployed  and  in  fact,  no  target  area  can  be  considered  to  be  inaccessible.  Thus, 
it  becomes  necessary  to  determine  a  defense  which  wilk  be  effective  against 
enemy  attacks  of  xinknown  duration,  strength,  location  and  time.  It  may  happen 
that  through  analysis,  subject  to  the  hypothesis  considered,  the  supply  of 
defensive  weapons  is  not  large  enough  or  the  defensive  weapon  system  has 
excessive  target  acquisition  times  to  defend  the  target  areas  adequately.  If 
so,  this  indicates  areas  of  weakness  that  should  be  remedied  either  thro'ugh 
increased  production  or  improved  engineering. 

In  the  following  pages  we  shall  describe  a  method  of  determining  how  the 
defensive  weapons  should  be  allocated  to  the  various  defense  areas.  At  the 
same  time  optimal  attack  policies  will  be  determined  and  the  best  defense 
against  these  optimal  attacks  will  be  determined.  In  the  assumption  involved 
we  shall  be  sometimes  optimistic  and  sometimes  pessimistic  about  the  capabili¬ 
ties  of  the  weapon  systems  of  both  sides,  the  intelligence  we  have  concerning 
the  location  of  targets,  their  worth,  the  number  of  the  defensive  weapons  and 
the  strength  of  the  attack. 

The  solution  to  the  problem  of  the  "best"  defense  will  be  somewhat  of  a 
compromise  solution  since  we  will  assume  that  the  enemy  wishes  to  do  his  best 
and  the  defense  will  attempt  to  do  his  best  against  the  attack.  Thus,  the 
defense  chosen  will  minimize  the  effect  of  the  attackers  best  offense  against 
any  assumed  defense.  This  defense  has  the  possible  drawback  effect  that  if 
the  attacker  does  not  use  an  optimal  attack  policy  (against  some  assiimed 
defense)  there  may  exist  other  defenses  which  will  result  in  less  expected 
damage.  We  shall  elaborate  on  this  point  after  formiilating  and  indicating 
the  method  of  solution  of  the  problem. 
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Problems  involving  optlmuin  deployment  of  attacking  missiles  against 
targets  which  are  defended  by  a  known  number  of  defensive  missiles  have  been 

1  o 

considered  by  Perkins'^  euad  Lindsey^".  While  the  problem  posed  here  is  somewhat 
related  to  the  problem  considered  by  Perkins  and  Lindsey,  the  computational 
procediore  and  localization  of  the  problem  differ.  W'e  shall  indicate  in  the 
following  how  the  problem  considered  by  Perkins  and  Lindsey  may  be  reformulated 
and  solved  using  the  methods  adopted  here. 

STATEMENT  OF  THE  PROBLEM 

Suppose  we  have  n  targets  which  are  to  be  defended  by  a  total  number,  D, 
of  defensive  weapons  (e.g.,  antimissile  missiles)  against  a  total  number,  A, 
of  offensive  weapons.  We  make  the  following  assumptions: 

(a)  The  targets  are  independent.  This  implies  that  once  the  de¬ 
fensive  weapons  are  allocated  to  specific  targets  they  cannot  be  reassigned 
to  other  targets  during  an  attack  (or  an  imminent  attack)  and  that  attacking 
weapons  that  miss  their  assigned  target  do  not  strike  another  target  in  the 
teurget  complex, 

(b)  The  target  locations  are  known  to  the  offense  and  each  side 
asBicms  the  seme  worth  to  a  particular  target. 

(c)  The  total  number  of  defensive  and  offensive  weapons  are  known 
to  each  side.  Their  distribution  among  the  targets  is  not  known. 

(d)  The  defensive  forces  assume  that  the  offense  will  deploy  their 
weapoai  in  such  a  way  as  to  maximize  the  expected  damage.  That  is,  if  the 
offense  knew  of  the  distribution  of  the  defensive  weapons  among  the  target 
areas  then  it  assigns  the  attacking  missiles  in  such  a  way  as  to  maximize  the 
expected  damage. 

(e)  If  a  missile  eludes  the  defense  then  the  probability  of  de¬ 
struction  of  the  target  Is  1. 

(f)  The  probability  that  the  defensive  weapons  fall  to  destroy  one 
or  more  of  the  attacking  missiles  (against  the  ith  target)  Is  a  known  function, 
p,  of  the  number  of  attacking  missiles,  the  number  of  defensive  missiles 
assigned  to  the  1th  target  and  the  probability  that  a  defensive  weapon 
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intercepts  an  attacking  missile  (if  assigned  to  it).  The  probability  that  a 
defensive  weapon  intercepts  an  attacking  missile  is  a  constant  k,  the  same 
for  all  defensive  weapons.  Fiarther  p  =  p(d^,a^>k)  has  the  properties  that 
p(0,a^,k)  =  1  if  a^  >  0  and  p(0,0,k)  =  0. 


(g)  Let  w^  be  worth  of  target  i;  then  the  expected  damage  to  the 
n 
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(h)  The  defense  philosophy  is  to  deploy  the  defensive  weapons  in 
such  a  fashion  so  as  to  minimize  the  expected  damage  resulting  from  the 
attackers  optimal  policies. 


The  problem,  then,  is  to  choose  a  set  of  <3-^'^  so  to  minimize 


n 


Max 


jS  }  .  ^ 

n  1=1 


^  w.p(d.,a.,k)  ,  i.e.,  we  wish  to  determine 


Min 

{d^ ,  . .  . ,  d^  } 


r  Max 


il-.  4-/^  r\  rs-r\f\  n  +  n  /^t*i  e  +  V»o  + 


w^p(d^,a^,k) 


^  1  n  1=1 


n  n 

\  —  7*^ 
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DYNAMIC  EROGRAMMING  FORMULATION  OF  THE  PROBLEM 


Our  goal  is  to  provide  a  feasible  computational  procedure  to  evaluate 

G(d,  ,...,d  ;a^,...,a  )  =  Min 
'1  ^  n  1'  n  ,, 


n 


(di,...,dn} 


n 


Max  )  w  p(d  ,a  ,k) 

{a  ,  ...,a  }  /-d 
.1  n  1=1 

n 


(2) 


subject  to  the  conditions  that  A  = 


a.  ,  D  = 
1 


d.  where  A  and  D  are 
1 


i=l 


i=l 


fixed  non- negative  integers. 
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To  obtain  the  computational  procedure  we  consider  A,  D  and  n  to  be 

parameters.  In  this  way  we  shall  be  able  to  reduce  this  problem  to  a  sequence 

of  simpler  problems  or  stages  and  at  the  ith  stage  we  will  be  able  to  determine 

a  specific  allocation  d^  to  the  ith  target.  To  attain  this  simplification  we 

use  the  functional  equation  technique  of  dynamic  programming  and  Bellman's 

3 

Optimality  Principle  .  First  we  note  that  the  value  of  G  will  depend  on  A, 

D  and  n  hence,  we  let 


N 


fj^(D,A)  =  Min 


or 


f^(D,A)  =  Min 


{d^, . . . ,d^} 


Max 


Max 


)  w.p(d.,a.,k) 

-Vi4  "  "  " 


for  1  <  K  <  n  (3) 


N-1 


{a^,  •  •  •  ^  Spj ) 


i=l 

1  <  N  <  n 


(^) 


Now  (4)  can  be  written  as 


f  (D,A)  =  Min 


N 


(5) 


To  obtain  (5)  from  (4)  we  have  employed  Bellman's  Optimality  Principle'^. 
For  N  =  1  we  have 


f^(D,A)  =  w^p(D,A,k)  (6) 

which  implies  that  for  a  constant  A,  p  decreases  as  D  increases  and  for 
constant  D,  p  increases  with  increasing  A.  Stated  differently  it  means  that 
if  we  only  have  one  target  to  defend  (or  attack)  then  we  defend  (attack)  it 
with  our  total  resources. 

Examination  of  (5)j  particular  the  quantity  within  the  braces  shows 

that  as  a  preliminary  calculation  we  are  to  calculate  optimal  attack  policies 

for  every  defense  of  the  Nth  target.  This  is,  in  essence,  the  problem 

2  2 

considered  by  Perkins  and  Lindsey  .  For,  if  we  drop  the  Min  operation  in 
equation  (2)  and  assume  that  D  is  fixed  and  each  is  known  the  maximum 
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value  of 


r\ 

WiP((ii,a^,k)  depends  only  on  A  and  n.  Hence  defining  6jj(A) 
i=l  n 

as  the  max  value  of  )  w  p(d.,a.,k)  we  have 

L _ I  1  11 

i=l 


Max  {  w^p(a. 
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/  A  ^  \ 


anfl 


_/a\  1.-\ 
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The  approach,  indicated  by  equation  (7)^  bas  the  advantage  of  compu¬ 
tational  simplicity  and  allows  one  to  assess  the  effect,  on  the  optimal  attack 
policy,  resulting  from  the  addition  of  more  targets  and  the  change  in  A. 


r'OMPTrPA'PTnMAT  PRrr'ITTITTRm 

V^X.iC~l.a.  U-l  »  I  t 


The  dynamic  programming  formulation  imbeds  the  problem  within  a  family 
of  analogous  problems  in  which  the  parameters,  D,  A  and  n  assme  sets  of 
values.  Examination  of  equations (5)  and  (6)  shows  that  the  first  step  in  the 
computation  is  to  assign  a  grid  of  values  for  D,  A  and  n.  Since  D,  A  and  n 
are  non- negative  integers  then  an  obvious  choice  wo\ild  be  to  let  D,  A  and  n 
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A  careful  analysis  of  equations  (5)  and  (6)  showsthat  in  the  process  of  calcu¬ 
lation  we  will  obtain  solutions  to  a  variety  of  sub  problems,  that  is,  problems 


wiiich.  msy  involvs  sius.Hoi’  nunib0i*s  of  ‘tai'gc'bs^  dsfonsivc  and  offsnsivc  woapons^ 
Furthermore  we  use  the  solutions  for  these  problems  to  obtain  solutions  of 
problems  of  greater  complexity. 


The  choice  or  selection  of  grid  values  for  A  and  D  will,  in  general, 
depend  on  the  magnitudes  of  A  and  D  and  also  the  memory  capacity  of  the 
available  digital  computer  as  well  as  the  time  and  cost  requirements  to  obtain 


OXXC  DL^J-U  UJ.v^xl^  ^  iiu  C2.j_  x  u.a.^0  v^cui  x  wj. iiiL4u.c* 


iH  mi  T“  Ci  V^TV\  *1  Q 

VX^'  UV...X  ilLXXl^ 


-LU. 


size  should  be  selected. 


wiiue  uxic  (jxiuxce  ui  ^xxu.  v<a.x.uc:c>  xux'  u  cuiu.  xi  iicid  uccii  iiidu-c  uxic:  'jctx’J uxci uxvjxi 

starts  with  equation  (6).  That  is,  if  &  and  a  are  arbitrary  elements  of  the 
D  and  A  grid  then  f^(6,a)  is  computed,  stored  and  printed  for  each  6  and  a 
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satisfying  the  relations  0<6<D,  0<a<A.  The  set  of  values,  {f^(6,a)), 
is  then  used  to  compute  f2(5>Q:)  from  equation  (5)-  For  N  =  2  equation  (5) 
becomes 


f  (6,0:)  =  Min  {Max 
d2 


W2p(d2,ag,k)  + 


■iCS-d^jO-a^) 


The  expression  h(d2)  =  Max 

^2 


+  f^(6-d2,  O-a^) 


is  a  function  of 


d2  and  we  must  determine  the  d^,  0  <  d^  <  6,  which  minimizes  h(d2).  The  value 
d^  is  a  function  of  6  £ind  o  and  this  value  gives  the  allocation  of  defensive 
weapons.  The  allocation  of  defensive  weapons  to  target  1,  is,  of  course, 
d^  =  b-dg.  The  set  of  values  {12(6,0)  j  0  <  6  <  D,  0  <  o  <  A)  is  stored  over 
the  set  of  values  of  fj^(6,a)  which  are  no  longer  needed  in  the  calculation. 
Also,  the  sets  of  values  12(6,0),  d2(6,a)  and  32(6,0)  are  printed  (and  stored 
if  Gonv'enient ) . 


The  above  procedure  is  repeated  until  f^(6,o)  is  obtained.  If  the  only 
case  of  interest  for  the  n  targets  is  when  6  =  D  and  o  =  A  then  this  calcu¬ 
lation  does  not  have  to  be  repeated  for  each  grid  point. 


The  structure  of  the  process  allows  the  decision  maker  to  assess  the 
effect  on  the  expected  damage  due  to  changes  in  the  number  of  targets,  de¬ 
fensive  and  offensive  weapons,  hence  the  solution  may  be  used  as  an  aid  for 
planning  purposes.  By  this  method  of  calculation  we  have  obtained  solutions 
for  1,  2  or  n  targets  for  all  ordered  pairs,  (&,a)  where  6  and  a  are  non¬ 
negative  integers  such  that  0  <  6<  D,  0  <  a  <  A.  Hence,  if  at  the  present 
time  the  defense  has  defensive  weapons,  the  offense  has  A^  offensive  weapons 
whereas  in  m  months  the  expected  numbers  of  defensive  weapons  will  be  and 
the  expected  number  of  offensive  weapons  will  be  A^  then  the  solution  gives 
the  expected  damage  for  each  case. 

Further  generalization  (or  restriction)  of  the  problem  is  possible  by 
modification  of  the  assumptions  (a)-(h).  For  example,  assumption  (e)  is  not 
necessary  or  realistic  for  most  conventional  weapons;  that  is,  it  may  not  be 
possible  to  totally  destroy  an  undefended  target  with  the  available  supply  of 
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offensive  missiles.  In  an  earlier  paper  ,  the  authors  considered  optimal 
attack  policies  against  targets  for  which  the  total  destruction  was  not 
generally  possible.  Each  target  had  a  survival  probability  associated  with 
it  which  could  be  interpreted  as  a  ftinction  of  its  defense  and  the  particular 
weapon  system  directed  against  it.  In  this  earlier  paper  the  resources  or 
supplies  of  the  offense  were  limited  and  included  several  classes  of  offensive 
vreapon  systems. 

The  problem  of  decoys  has  not  been  mentioned  above  but  may  be  included 
by  the  introduction  of  proper  constraints  or  by  redefining  the  function 
p(di^ai,k).  If  the  defensive  weapon  system  has  rapid  response  and  short  target 
discrimination  times  decoys  may  pose  no  problem;  however  since  cases  for  which 
this  is  not  true  exist  any  realistic  problem  should  include  a  treatment  of 
decoys.  Perkins^  includes  these  features  in  his  discussion. 


In  the  following  pages  we  consider  a  simple  numerical  example  which 
illustrates  the  method  of  solution  and  information  obtained  by  solving  equations 
(5)  and  (6)  by  the  above  described  method.  The  straight  forward  solution  of 
equations  (5)  and  (6)  as  described  above  requires  the  storage  of  a  table  of 
(D+1)(A+1)  elements  in  order  to  proceed  to  the  next  stage  of  the  calculations; 
hence  the  magnitude  of  D  and  A  which  can  be  handled  is  a  function  of  the 
available  high  speed  memory  of  the  available  computer.  In  our  example  this 
poses  no  problem  since  we  are  interested  not  in  solving  the  most  general  problem 
but  in  illustrating  the  method  of  calculation  and  the  kind  of  information  ob¬ 
tained  from  this  method  of  solution.  In  our  example  p  =  p(d.ja.,k)  is  the 

2  ^  ^ 
function  used  by  Lindsey  ,  however,  d^  is  not  known,  a  priori,  but  must  be 

determined.  As  stated  previously,  this  determines  the  "best"  defense  against 

the  attackers'  optimal  attack  policies, 

AN  ILLUSTRATIVE  EXAMPLE 


Let  p  =  p(d.,a.,k)  =  1  -  i  1-exp  (  -kd.  /s'"  } 

=  0 

A  =  10,  D  =  25,  k  =  1/2 


a . 

1 


if  a.  >  0 
1 

if  a.  =  0 
1 


and  suppose  we  have  three  targets 

w,  =  16,  w  =  L 


J  )i  -TT-  _  '-7 

Lh-  a-iiu.  —  (  . 


T^,  T^  whose  respective  worths  are 
Our  goal  is  to  evaluate  the  function  given  by 
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equation  2.  To  accomplish  this  we  imbed  the  problem  into  a  family  of  analo¬ 
gous  problems,  that  is,  we  consider  the  numbers  D  and  A  to  be  parameters  and 
allow  the  niomber  of  targets  to  vary  also.  We  let  6  and  a  be  any  arbitrary 
non- negative  integer  less  than  or  equal  to  D  and  A,  respectively  and  start 
the  calculation  with  the  computation  of  f^(6,a)  (equation  6)  for  all  admissi¬ 
ble  values  of  &  and  a.  This  calculation  gives  us  the  expected  damage  of  target 
T^  when  the  defense  consists  of  &  missiles  and  the  attacker  has  a  missiles. 

The  results  of  these  computations  are  stored  in  the  high  speed  memory  of  the 
computing  device  and  are  given  also  in  Table  I.  Since,  for  one  target,  we 
defend  with  the  total  number  of  defensive  weapons  available  and  the  attacker 
attacks  with  what  he  has  available  there  is  no  real  need  to  print  the  corre¬ 
sponding  defense  and  attack  policies.  The  next  step  is  to  compute  f„{6,a)  for 
all  admissible  6, a  from  equation  (5),  with  N  =  2,  l.e.. 


f  (6,0:)  =  Min 


Max{w2p(d2,ag,k) 

a„ 


+  f3_(8- 


d„,  a-a_)} 


(5) 


The  calculation  of  f2(6,o;)  starts  with  6  =  D  and  a  =  A  and  the  result  f2(D,A) 
is  stored  over  f^(D,A)  which  is  no  longer  needed  in  the  calculation  process. 
The  values  f2(D,A),  1^(0, A)  and  3^(0, A)  are  also  printed.  Then  f2(D-l,A)  is 
computed  and  stored  over  f^(D-l,A)  and  f2(D-l,A),  d2(D-l,A),  a2(D-l,A)  are 
printed.  The  process  is  continued  until  the  complete  set  of  values  f2(6,Q:) 
is  computed.  At  the  completion  of  this  stage  we  have  obtained  solutions  for 
cases  involving  2  targets,  5  defensive  weapons  and  a  offensive  weapons.  The 
expected  damage,  the  allocation  of  defensive  weapons  and  the  optimal  attack 
policy  are  presented  in  Tables  II,  III  and  TV.  Thus,  for  example,  for  2 
targets,  20  defensive  weapons  and  6  attack  weapons  we  find  in  Table  II  that 
the  expected  damage  is  lh.66.  Similarly  we  find,  in  Table  III  that  target  T^ 
should  be  defended  with  7  defensive  weapons,  and  target  T^  should  be  defended 
with  20-7=13  defensive  weapons.  The  optimal  attack  policy  which  yields  this 
of  expected  damage  is  given  in  Table  IV,  i.e.,  a^  =  0  and  a  =  6  -  =  6. 

The  next  stage  of  the  computation  is  to  compute  fj(6,a)  from  equation  (5) 


It 


with  N=3,  i.e. , 


Min 

d. 


Max  {w^p(d^,a^,k)  +  f2(&-d^,  a-a^)) 

^3 


The  results  of  this  calculation  are  presented  in  Tables  V,  VI  and  VII.  At 
this  stage  the  above  example  has  been  solved  as  well  as  many  other  similar 


pX  WUJ_c:iuo  f 


^  r\  1  ^>nc< 


The  entry,  in  Table  V,  corresponding  to  D  =  23,  A  =  10  is  expected  damage 
on  three  targets  which  are  defended  by  25  defensive  weapons  against  10  attack¬ 
ing  weapons.  We  find  that  f^(25,  lO)  =  25.35.  The  defense  of  T^  and  the 
optimal  attack  policy  against  T^  are  given  in  Tables  VI  and  VII  respectively. 
Thus,  we  find  that  d^  =  2  and  a^  =  2  which  leaves  23  defensive  weapons  and  8 
attack  weapons  to  be  allocated  to  T^  and  T^.  From  Table  III  we  find  that 
*^2  ^  ^  from  Table  IV  a^  =  a^(2^,d)  =  4  hence  d^  =  l4  and  a^  =  4. 

Thus,  the  best  defense  for  the  3  targets  (with  their  25  defense  weapons) 
against  the  attackers  10  offensive  weapons  is  (d^,d2,d^)  =  (14,9^2).  This 
defense  of  the  three  targets  Insures  that  the  expected  damage  will  not  exceed 
25.55,  The  attack  policy  which  yields  this  level  of  expected  damage  is 
^^l’^2’^5^  =  (4,4,2)  and  any  other  attack  policy  against  these  three  targets 
and  the  given  defense  will  result  in  a  lower  value  of  the  expected  damage.  As 
a  further  illustration  of  the  above  suppose  there  exist  only  20  weapons  to 
defend  the  three  targets  and  the  attacker  has  six  offensive  weapons .  Then  we 
find  that  f^(20,6)  =  17.43  and  the  best  defense  is  (10,7,3).  The  optimal 
attack  policy  then  is  (0,4,2).  If  the  attacker  allocates  his  missiles  such 
that  (a^,a2,a^)  =  (3,2,1)  against  these  three  targets  with  the  defense  (lO,7,3) 
then  the  expected  damage  will  be  13.47.  Furthermore,  if  for  some  reason,  the 
above  defense  is  not  acceptable  to  the  defenders,  and  should  be,  say  (9^8,3) 
then  this  results  in  a  penalty  for  then  the  maximum  expected  damage  will  be 
greater  than  17.43,  viz.,  17*72  for  an  attack  policy  of  (4,0,2)..  This  infor¬ 
mation  is  not  generally  printed  out  but  is  obtained  actually  in  the  course  of 
the  calculation.  This  information  may  be  useful  in  some  instances  when  we  have 
associated  cost  functions  for  then  it  may  be  advantageous  to  have  near-optimal 
policies  of  defense.  The  method  of  solution  gives  results  which  may  be  useful 
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for  many  purposes,  for  example,  one  may  construct  iso-damage  curves  from  the 

f  -tables  and  thus  be  able  to  anticipate  future  needs  of  defensive  weapons  as 
n 

the  attachers  supply  of  offensive  missiles  increases.  Furthermore,  from  the 
f^-tables  one  can  estimate  the  effect  on  the  expected  damage  which  results 
from  over  or  under- estimating  the  number  of  attack  missiles  available  to  the 
attacker. 

It  may  happen,  that  because  of  transportation,  cost,  defense  site  pleui- 
ning  or  structure  of  the  defense  system  one  does  not  desire  the  allocation  of 
an  individual  missile  to  a  target  but  rather  wishes  to  allocate  the  defensive 
weapons  in  blocks  of  m.  If  this  is  desired  then  the  grid  of  D-  values  should 
be  chosen  to  be  0,  m,  2m,  . . . ,  D  ■  Jm.  In  fact,  such  a  method  actually  simpli¬ 
fies  the  calculation,  reduces  memory  recjuirements  and  computation  time.  Further¬ 
more,  the  method  of  solution  described  above  allows  one  to  assess  the  effect 
that  such  an  allocation  has  on  the  expected  deunage  and  yields  Information  as  to 
what  the  size  of  the  allocation  block  0houj.d  be,  That  is,  one  actually  performs 
a  sensitivity  analysis  of  the  problem  in  the  course  of  solution. 

W.  J.  SAGGO 

R.  E.  SHEAS 
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